Full versus first stage replica symmetry breaking in spin glasses. 
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A short survey is presented on spin-glass-like states characteristics in complex nonmagnetic sys- 
tems. We discuss the interplay of the interaction structure and symmetry with the classification 
scenarios of the replica symmetry breaking. It is shown that the kind of the transition to the non- 
ergodic state depends not only on the presence or absence of the reflection symmetry but on the 
number of interacting operators and their individual characteristics. 
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I. INTRODUCTION 

The reflection symmetry plays the crucial role defining 
the character of phase transition in nonrandom mean- 
field (MF) models.EJ Generally speaking the presence of 
the terms without reflection symmetry usually results 
in the first order phase transition, while in the absence 
of such terms the transition is of the second order. In 
the case of random MF systems the absence of refiection 
symmetry also leads to a special form of the free energy 
functional that differs from the symmetrical case. As a 
consequence, the scenarios of replica symmetry breaking 
(RSB) are different for these two cases. However, not 
only symmetry determines the transition to nonergodic 
state. Extending the class of models permits considering 
the role of different factors in the scenarios of appearance 
of SG-type states. In this paper we try to use our recent 
results for different models for spin-glass-like states in 
complex nonmagnetic systems (see Ref.g for a review) 
to investigate how the interaction type correlates with 
the spin glass (SG) behavior. 

The theory of spin glasses developed as an attempt to 
describe unordered equilibrium freezing of spins in actual 
dilute magnetic systems with disorder and frustration. 
This problem was sooru solved at the mean- field level 
Edwards and AndersorJp Sherringt(m,and Kirkpatrick,Q 
Almeida and ThoulessJ3 and ParisiErEl [see Ref. ^ for a 
review]. The Sherrington-Kirkpatrick approach to the 
spin-glass theory starts from the Hamiltonian 



if = -i^J„- [/.{/, 



(1) 



i^J 



It describes Ising spins U located on the lattice sites i. 
The quenched interactions Jij are distributed with the 
Gaussian probability. 
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where J = J/y/N , Jo — Jq/N and N s the number 
of sites. To perform averaging over disorder in this case 
one has to average the quenched free energy F rather 
than the partition sum Z itself. Such averaging is usually 



performed using the replica method. After that the free 
energy becomes the function of the order parameters that 
depend on replica indices: 



F^F{x°',q°"'), 



N N 



(3) 

(4) 



The free energy F{x" ,q"") has an extremum for the 
replica symmetric (RS) solution when all q"" are equal. 
However this state is unstable under RSB. Parisi pro- 
posed the method of RSB step by step with the limit full 
RSB (FRSB) when q"^ becomes a continuous function 
of a parameter x. This approach allows to describe the 
main features of the experiments on spin glasses. Namely, 
in the framework of the equilibrium approach, the spin 
glass phase with qualitatively correct boundaries was ob- 
tained and the difference in the behavior of magnetic sus- 
ceptibility in field-cooled and zero-field-cooled cases was 
explained. 

So, the problem of theoretical description of SG per 
se was solved in principle and during that time different 
other models appeared without any connection to real 
experiments and real physical systems. The main fea- 
ture of these models was the absence of time reversal 
symmetry - in contrast to the SK model. The most in- 
vestigated model among those are the p-spin models and 
the Potts models, considered, for example, in Refs. ^ 
p^ . The spherical p-spin model [see, e.g., Ref. |l^ was 
believed for a long time to be a generic for this class 
of models. From the point of view of RSB the main 
feature of this model is the stability of the first step of 
RSB (IRSB) down to zero temperature. Also, the or- 
der parameter behaves stepwise. Although this model 
was not aimed to describe any actual glass it appears 
to be very interesting because its behavior gives a sce- 
nario for real liquid-glass transition: two critical temper- 
atures, the number of metastable states similar to that 
obtained in numerical modeling. It should be noted that 
the structure of the dynamical equations for the corre- 
lation functions are identical for the supercooled liquids 
in the mode^coupling theory and for the p-spin spherical 
SG modelta 



Based mainly on the investigations of these two models 
- SK and p-spin spherical - a conclusion appears in the 
literature attributing the two classes of universality to 
the models with and without reflection symmetry. In 
the disordered case a number of attempts were made to 
formulate a kind of universality rules based on the mean- 
field investigatipn of the model systems with the random 
interactions. BE3 

Looking now in general at the free energy series over 
the glass order parameter we see that the series contain 
explicitly the terms which can be classified by the reflec- 
tion symmetry, 



AF^ 



lini i y L. + a-i5q°''^5qf^^5q'°'+ 
n-s-O n ^ — ' 



NkT n^u n 

a'^iSq°'^Y + a4Sq"^Sq'^''6q'>^Sq^°'...] , (5) 

and the part without the reflection symmetry: the terms 
with three identical replica indices: 



NkT 



\im-y[... + b3iSq''^f + ...+ 

bi6q°'^Sq'^'^Sq-i°'Sq^°'...] . (6) 



Thus, a natural question arises: whether there can be 
made a general statement regarding the behavior of SG 
models with and without reflection symmetry? And do 
all models of the first type behave in fact as SK model and 
all models of the second type as p-spin spherical model? 
In this paper we try to answer this question. We prove 
that for arbitrary models with the reflection symmetry 
the Parisi FRSB always takes place. In the absence of the 
reflection symmetry the situation is not so definite and 
the behavior of the system depends on some additional 
characteristics. In any case it is not always similar to that 
of he p-spin spherical model, as it was usually believed: 
we give the counterexamples. 



II. GENERALIZED SK MODEL: FRSB 
A. SK-model w^ith the reflection symmetry 

In this case it occurs to be possible to prove a kind of 
a theorem. 

First, we consider a generalized model with reflection 
symmetry with the Hamiltonian (nh, with the interac- 
tions distributed according to Eq.O) and with the ar- 
bitrary diagonal operators U. The reflection symmetry 
implies that for any integer k, 



Tr 



(7( 



2fc+l) 



0. 



(7) 



The saddle point conditions for the free energy aver- 
aged over disorder produces the glass order parameter 



g"'^ ^ Tr [U°'U'^ exp (9)] /Tr [exp (6)], 



and the auxiliary order parameter 

w°' = Tr [(C/")2 exp (6*)] /Tr [exp (9)]. (9) 

Here 

= ^^w"(L/")2+t2^g"^[/°C/'3, (10) 



a>P 



where t = J/kT and we choose Jo = for simplicity. 

In the RS-approximation we find the (trivial) solution 
9rs = 0. The bifurcation condition looks like in this case: 



l-tlw^{Q = 0. 



(11) 



This equation coincides with Arepi(RS) = [see, e.g., 
Ref. @|. It is very important that it is zero solution that 
bifurcates. 

Investigating IRSB, 2RSB, 3RSB, and so on, we see 
that the equations for the glass order parameters always 
contain the quantity 



Tr[C/exp(0„RSB)]/Tr[exp(0„RSB)]. 



(12) 



Therefore, one of the solutions of these equation is triv- 
ial at each of the RSB-steps, and the appearance of the 
7iRSB solution can be regarded as the bifurcation of the 
trivial (n — 1)RSB solution. In this case, the equation 
AnRSB = coincides with the corresponding branching 
condition (pi]). This means that in any case, the nRSB 
solutions at different stages of the symmetry breaking can 
exist at the temperature T < Tc determined by this bi- 
furcation condition, and so we always can look for FRSB 
solution. Writing the free energy as a series over Sq"/^ 
near Tc (up to the fourth order of magnitude inclusively) 
we obtain q{x) = ex in the leading approximation [sim- 
ilar procedure in details was described in |lj] . It is also 
possible to write the free energy in the form of Parisi 
with the only difference in the boundary conditions for 
the Parisi function ch that now reads: 



(l,y) =lnTr<^ exp 



tyU+*-{w-q{l))U^ 



(13) 



(8) 



Thus we have shown that in the case of systems with 
reversal symmetry, the infinite FRSB occurs at the very 
point at which the RS solution becomes unstable. In par- 
ticular, our result means that magnetic systems of arbi- 
trary spin with the interaction between the z-components 
behave in the same way. 



B. SK-model without the reflection symmetry 

We consider below the models without the reflection 
symmetry. [These models correspond to some real phys- 
ical systems.] It implies that we have the Hamiltonian 
(E|)-(p|) but without the condition (R) for the operators U. 
It is easy to trace how the proof given above fails using 
the model proposed in Ref. W^ 



The difference between two cases is already manifested 
in the RS approximation. In the case when the condition 
(0) is not fulfilled for the Hamiltonian (|]) there is no 
trivial solution for the order parameters. The disorder 
smears out the first-order phase transition; hence, in- 
stead of a transition, there is a smooth increase in the 
order parameters (both glass and regular) as the temper- 
ature decreases. This situation is seen in experiments on 
orientational glass phaseJp ortho— para-hydrogen mixed 
crystals and in Ar — N2 .113 These substances present mix- 
tures of spherically symmetric molecules and momentum 
bearing molecules. The corresponding glass was investi- 
gated on the base of the Hamiltonian (^) with U = Q, 
where Q = 3 J^ - 2, J = l.llZl The RSB so ution branches 
continuously and smoothly on cooling breaking the RS 
results in a transition to the nonergodic phase of the 
quadrupolar glass. 

Another example of a SG-like phase in the molecu- 
lar crystal i^he pure para — H2 (or ortho — D2) un- 
der pressure.tS The possibility of the orientational order 
in the systems with the initially spherically symmetric 
molecule states is due to the involving of higher order or- 
bital moments J — 2,4... under pressure. With increase 
of the density the anisotropic interaction potential and 
the crystal field grow rapidly and the energy of the many 
body system can be lowered taking the advantage of the 
anisotropic interactions. The long range orientational or- 
der appears abruptly at a fixed value of pressure through 
the first order phase transition just as it takes place in 
ortho-para mixtures when the concentration of moment 
bearing molecules achieves certain fixed value. In the in- 
termediate concentration range the frustration and dis- 
order motivate the investigation of the quadrupole glass 
with J — 2. Such a theory was constructed in |l^. The 
essential feature of the obtained intermediate phase is the 
coexistence of the orientational glass phase with the long 
range orientational order as it is seen in the experiment. 

We consider two other models describing SG-like 
states in real complex nonmagnetic systems, namely, in 
the cluster systems. Although they are not the mixtures 
of the different kinds of particles with different interac- 
tions, one can find frustration and disorder, that is the 
background to consider these systems in the spirit of SG 
theory. Now the operator U in (nf) is replaced with con- 
tinuous functions of the angles. 

In Ref. ES a model for the low-temperature transition 
to the orientational glass state in solid molecular Cgo was 
developed. Although the molecules have nearly spheri- 
cal shape, at low temperature there are two pronounced 
minima in the anisotropic part of the intermolecular in- 
teraction energy. It is possible to trace an analogy with 
the mixtures studying the role of the different types of 
the mutual molecular orientations. As a result, a model 
is constructed where the role of spin is played by certain 
combinations of the cubic harmonics. The results agree 
well with the experimental data: the coexistence of the 
glass state and the long-range orientational order and the 
existence of a wide maximum on the curve for the orien- 



tational part of the heat capacity. Moreover, the above 
model permits considering tlie pressure dependence of 
the orientational transitions. Ell 

The other model we would like to mention is the SG- 
like freezing of clusters of different symmetries in super- 
cooled liquids that gives a possible description of liquid- 
glass transition. In Ref.E2 we use the microscopic ap- 
proach based on the equations for the distribution func- 
tions which in spirit of Bogoliubov hierarchy give us an 
opportunity to analyze the intercluster interaction. We 
show that there exists a region of densities and temper- 
atures where this interaction changes sign as a function 
of the cluster radius and there is hence frustration in the 
system. This is the base to write a Hamiltonian of the 
form (Q) with different point group harmonics for U and 
use standard methods of SG theory to describe the real 
glasses. 

So, we have considered a set of models with two- 
particle interaction where the absence of the refiection 
symmetry is caused by the characteristics of the opera- 
tors U themselves. In this case the RSB-solution bifur- 
cates from the RS -solution smoothly, without a jump, 
and the coexistence of the glass order with the long range 
regular order takes place. 



III. GENERALIZED p-SPIN MODEL 

We consider now a generalization of the well-known Jo- 
spin modelll3 of Ising spins where spins are replaced by 
arbitrary diagonal operators. Then the Hamiltonian is 



H = - 



E 

ii<i2...<ip 



^i\...ip^ i\ ^ i2 •••^ij 



(14) 



where i — 1,2, ...N, p is the number of interacting par- 
ticles and U is an arbitrary diagonal operator such that 
TrC/ = 0. We do not specify its form here, in order to 
use general expressions further. The independent inter- 
actions have the Gaussian distribution 



P{Jii...ip) = — ^-^~ exp 



p!J2 



(15) 



Using the standard procedure of replica approach we 
obtain the free energy and the equations for the order 
parameters (ph and (0), but now with. 



t^ 



0=p^^(g"/3)(p-l)c/O[//3 + 



a>l3 



PTE(^")^'"'^(f^")'- (16) 



We perform the first stage RSB [n replicas are divided 
into n/mi groups each containing mi replicas] and obtain 




FIG. 1. (Color online) Temperature dependence of the glass 
order parameters for the quadrupole glass with three particle 
interaction for J = 1. The RSB occurs at the temperature 
corresponding to the condition Aj^sj^cpi = 0. 



the free energy in the form 
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FIG. 2. (Color online) Temperature dependence of the order 
parameters for the quadrupole glass with three particle inter- 
action for J — 2. The transition RS — IRSB takes place at 
the point defined by the condition m — 1. The glass order 
parameter vi has the jump at this point. Here x is the regular 
orientational order parameter. 



i^iRSB = -NkT[mit''{p - 1)^ 



-t^p-l)^ 



— fdz^ln /ds'^[Trexp0iRSB]'"1. (17) 
mi J J 



Here q"" = ri if a and /3 are from the different groups 
and q"^ = ri + vi otherwise, 



/pri(P~i) ^^ 

71RSB = Zt\l U 



p[in+vit'^-n(p-^)] ^ 



+ .^^[-^""(-^+-^)""V-. (18) 



2RSB appears]: 

A(iR,sB)rcpi == 1 - jPiP - l){ri + vi^P-^^ X 
/"dz° /"ds°[Trexp0iRSB]"^x 



/ Tr([/^expgiRSB 
I Trexp6'iRSB 



Tr(C/exp6'iRSB) 



Trexp6'iRSB 
/"ds'^(Trexp0iRSB)"^ 



(19) 



So, in the case of the three-particle interaction be- 
tween quadrupoles with J = 2 as well as with J — 1, 
the first stage RSB is stable only in the finite region of 
temperatures and not down to zero temperature [as was 
supposed in Ref. |l2[. This is the key result. As concerns 
the models with the multiple interactions, this property 
was investigated for the Potts model with three states 
in Ref. |2^ [see, also Ref. p3 and for more complicated 
models with two interactions in Refs. E5^-E8l 



We performed the detailed calculations for two mod- 
els with p — 5 and where U were represented by the 
quadrupolar moment with J = 1,2. The results of the 
calculations are illustrated in Fig. ^ ( J = 1) and Fig. || 
(J = 2). The stability of the IRSB solution against fur- 
ther RSB was checked in the standard wayoEj looking 
at the positive values of A(iRSB)rcpi [defined as the bifur- 
cation point where the non-zero order parameter V2 in 



IV. CONCLUSION 

We considered the behavior of complex spin-glass-like 
systems. The set of physical systems can be divided into 
two classes depending on whether the reflection symme- 
try is present or not. 

We have shown that in the systems with the reflection 
symmetry the infinite FRSB takes place at the very point 



where the RS solution becomes unstable. This behavior 
is well known for the SK spin model. In particular, our 
result means that magnetic systems of arbitrary spin with 
the interaction between the z spin-components behave in 
the same way. 

If there is no reflection symmetry then the situation 
is not so definite. The behavior of the system depends 
on additional characteristics. An important property of 
such systems is the absence of a trivial RS-solution. We 
have considered a set of models with the two-particle 
interaction where the absence of the reflection symmetry 
is related to the structure of the [/-operators. In this 
case the RSB solution bifurcates from the RS solution 
smoothly, without a jump. The jump appears in the 3- 
quadrupole glass model with J = 2. The coexistence of 
the glass order with the long range regular order takes 
place in all the cases. 

The properties of the models considered in our paper 
are not similar to the properties of the p-spin spherical 
model as follows from three counterexamples. We have 



shown that in these cases under certain additional con- 
ditions there exists a finite domain of stability for the 
IRSB order parameters. This was apparently first shown 
for simple nonspherical models in Refs.^] and |l^ This ef- 
fect was discovered for the Potts model with three states 
in Rcf. ^ earlier. The FRSB is attained as a result of 
several successive transitions taking place as the temper- 
ature decreases. 
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